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Abstract. The Coma cluster exhibits evidence for a high-energetic non-thermal particle population. At frequencies 
> 1GHz recent radio halo observations confirm a significant spectral steepening of the volume-integrated emission. 
We calculate the volume- averaged high-energy spectrum due to inverse Compton scattering off the CMB radiation 
field and non-thermal bremsstrahlung according to an exponential cutoff in the electron spectrum as deduced from 
the radio observations. The synchrotron radiation from secondary pairs, created from the decay of charged mesons 
produced in hadronic pp-interactions, is found to set significant constraints on the energy content of relativistic 
hadrons in Coma. This limits the maximum flux at high energies. Our findings support a low ratio of relativistic 
hadron to thermal energy density. Predictions for Coma's high energy emission are discussed in the light of current 
and expected abilities of upcoming 7-ray instruments. 
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Clusters of galaxies are conglomerates of a large number of galaxies that are gravitationally bound and confining a 
large fraction of the mass in the universe. One of the controversially discussed properties of clusters of galaxies are 
their non-thermal components which include cosmic rays as well as turbulence and non-regular magnetic fields. The 
non-thermal pressure has an important impact on the evolution of galaxy clusters. To achieve a better understand- 
ing of these components is mainly driven by the importance of the non-thermal pressure in the evolution of galaxy 
clusters. In this paper we discuss only cosmic rays as the easiest testable one of all non-thermal components. Several 
mechanisms have been proposed that lead to relativistic particles in the intracluster medium (ICM), as e.g. particle ac- 
celeration during their formation and evolution, in e.g. merger shocks (Takiz grwa~fe Naito 20001 IGabici fe Blasi 2 003b, 
|Berrington fe Dermer 2003| IMiniati et al. 2001a|) . accretion shocks IjColafrancesco fc Blasi 1998|l . intergalactic termi- 
nation shocks from the winds of the galaxies IjVolk et al. 19 96). or reacceleration of injected mildly relativistic parti- 



cles from powerful cluster members (Enfilin et al. (1997) I. Indeed, the detection of synchrotron radiation from clus- 



ter radio halos and relics signals the existence of relativistic electrons (e~) in the intracluster medium (ICM). The 
most prominent cluster which possesses a radio halo is the Coma cluster (Abell 1656), located at redshift z=0.0232 
IjStruble fc Rood 1991(1 (corresponding to ~ 90 Mpc for Ho=75 km s _1 Mpc -1 ). At frequencies below ~ 1 GHz observed 
volume-integrated fluxes are satisfactorily fitted by a pure power law. Observations at larger frequencies gave evidence 
that a significant spectral steepening of the integrated emission occurs in Coma's radio halo (jSchlickeiser et al. (1987)1, 



recently confirmed bv lThierbach et al. 20031 (TH03'). A hard X-ray (HXR) excess has been detected by the Rossi X-ray 
Timing Explorer (RXTE) ( |Rephaeli et al. 1999| | and BeppoSAX (|Fusco-Femiano et al. (2004) I. Coma is one of the 



few clusters where also an EUV excess is conclusively established ( |Bowyer fc Berghofer (1998)|ILieu et al. 19991 . The 
spread of the soft X-ray to EUV emission is, however, still debated l|Kaastra et al. 2003J) . and may possibly extend 
up to the turnaround radius. We restrict here our considerations to the size of the radio halo (~10'x30'; TH03). By 
comparing the radio synchrotron spectrum with this excess radiation, interpreted as Inverse Compton (IC) scattering 
off photons from the cosmic microwave backgr ound (CMB) by the same e~ population, volume-averaged magnetic 
fields of B = 0.1 . . . 0.3^G have be en deduced (|Fusco-Femiano et al. (2004) |Rephaeli et al. 19"99| ). Faraday rotation 



measurements gave B ~ 2 . . . lQ/iG l|Kim et al. 19901 IClarke et al. 2001JI . 

Many models for non-thermal radiation from the Coma cluster predict significant emission at > 100 MeV due 



to both, relativistic e and ions IjSarazin 1999al |Atoyan fc Yolk 2000| |Miniati (200"3]) IGabici fc Blasi 20040 . They 
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often assume a straight power law to ~ 10 6 " 7 MeV for the e~ population responsible for the dominating synchrotron 
component. Here we investigate the consequences of the decline in the e~ spectrum at ~ 10 4 MeV as deduced from recent 
radio observations of Coma C for the expected high energy flux. In Sect. 2 the IC and non-thermal bremsstrahlung 
spectrum from the steepening e~ distribution is calculated. In Sect. 3 we derive limits imposed by the broadband 
observations for the 7r°-decay 7-ray component including its secondary pair initiated radiation. Finally we discuss 
Coma's detectability with current/future 7-ray instruments such as INTEGRAL, GLAST-LAT and the new Imaging 
Air Cherenkov Telescopes (IACTS). 



2. Non-thermal electron spectrum and radiation 



The volume-integrated radio emission from the radio halo has been studied in detail by e.g. Schlickeiser et al. (1987) 
IKim et al. 19901 |Venturi et al. (1990)1 |Giovannini et al. (1993)] |Deiss et al. (1997)| TH03. Fig. □ shows the volume- 
integrated radio continuum spectrum of the diffuse radio halo source Coma C as published in TH03 with the best fit 
model. TH03 confirmed the findings o f | Schlickeiser et al. (1987) that among the three basic models for cluster halos 
(primary electron model: Jaffe (1977)[ |Rephaeli (1979)| secondary el ectron model: e.g. IDennison 19801 in-situ accel- 
eration model: Jaffe (1977)[ Roland (1981)| |Schlickeiser et al. (1987) I the in-situ acceleration model fits the observed 
exponential steepening of the synchrotron spectrum best. This model, though discussed critically by Petrosian (2001) 
considers shock wave and resonant diffusion acceleration out of a thermal pool of particles where radiation losses and 
particle escape has been taken into account. A secondary origin of the radio halo has been pushed forward by many 
authors (e.g. IDennison 19801 |Blasi fc Colafrancesco (1999)| |Dolag fc EnBlin (2000)[ |Atoyan fc Volk 20001 IBlasi 20011 
IMiniati et al. 2001 bfl . Recently, however, arguments are given which suggest that secondary pairs as the underlying 
particle population of the radio halo emission are problematic IjBrunetti 20031 IKuo et al. 20 04). Along these lines 



Brunetti et al. (2004) found that the observations of non-thermal radiation of galaxy clusters are only reproducable 



within the picture of particle acceleration through cluster merger generated Alven waves, if the fraction of relativistic 
hadrons in the ICM is small (5-10%). This hadron content is insufficient to reproduce the radio halo from secondary 
pairs (see below). Curved spectra are also possible at an energy where losses balance the acceleration rate if the ac- 
celeration time decreases more slowly than the loss time. In the following we therefore consider an exponential shape 
of the e _ spectral distribution, suitable to explain the volume-averaged synchrotron spectrum , irrespective of its 
formation mechanism. This rather phenomenological ansatz will not shed light onto the mechanisms responsible for 
the formation of the e _ distribution, however, leads to model-independent constraints for the high-energy component 
arising from this leptonic particle population. 

We fit the radio flux density with a power law synchrotron spectrum extended by an exponential cutoff: 



7s yn (i/) cx v {3 s)/2 exp(-Vz7^) 

with s — 4.6 and v s = 0.44 GHz (TH03). Synchrotron-self absorption will affect the synchrotron spectrum at low 
energies. For an estimated path length of ~280 kpc through the cluster (TH03) we found the turnover frequencies 
at ~ a few 100 kHz (see Fig. QJ. Below ~ 0.2 MHz free-free absorption in the disk of the Milky Way suppresses the 
radio intensity from Coma observed at Earth. In addition, the Razin-Tsytovich effect causes the radio spectrum to 
decline rapidly to low frequencies from vr ^ 20(n e / ' Bq) Hz (n e is the e~ density in cm -3 , and Bq is the magnetic 
field in Gauss) with turning point between 0.01 . . . 0.6 MHz for the expected field strengths in Coma. Turnovers below 
MHz-frequencies are, however, not possible to detect with ground based radio observatories due to ionospheric effects, 
and have to await future space-based low-frequency observatories, for example a radio observatory on the Moon. 
For a given magnetic field B the corresponding volume-integrated e _ spectrum 

nip) = QoP~ s exp (-p/p c ) 



I Schlickeiser et al. (1987) ) with normalization Qq and cutoff e~ momentum p c can then be determined. Values for p c 
considering magnetic field strengths of 0.1. . . 6{iG lie at 2153 . . . 16667MeV/c. IC scattering off CMB photons with 
photon energies £cmb — 6 • 10 -4 eV (Tcmb = 2.7 K) by these e~ is therefore restricted to the Thomson regime. In 
the ^-approximation the intensity of the IC scattered radiation as a function of photon energy E 1 can be analytically 
expressed if the target photon density distribution n(e) is sufficiently peaked so that E c = (p c /m e c) 2 e rj (jp c /m e c) 2 Z: 
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with r the Gamma function, £ is Riemann's zeta function, m e c 2 the e rest mass, <jt = 6.65 



dr Coma's 



luminosity distance and n(e) is the CMB photon density. Fig. shows the resulting IC spectra for an equiparti- 
tion magnetic field ranging between 0.68 . . . 1.9/xG (TH03), for the central magnetic field in Coma C (B ~ 6/iG: 
IFeretti et al. 1995|l and for B = 0.1/xG appropriate to explain the HXR excess emission. 

The non-thermal volume-averaged bremsstrahlung intensity using the primordial 4 He mass fraction of 0.24 



1.18n. t c£ 7 



max(£ 7 m e c,p m i n ) 



dpp 2r n(p) 



da 
~dK 



is calculated in the relativistic limit using the differential cross section from Blumenthal & Gould (1970) For a mean 
gas density ni ~ 10~ 3 cm~ 3 in Coma we find the compound IC and bremsstrahlung spectrum from the observed 
synchrotron-emitting e - distribution always below the EGRET upper limit (Rcimc r" et al. 2003jl . The steepening of 
the e~ spectrum at 10 3 " 4 MeV causes the IC component to decline at ~1-10 MeV, non-thermal bremsstrahlung 
dominates till its decline at a few GcV. This is in contrast to works where the primary e~ spectrum extends to several 
10 7 MeV (e.g. |Atoyan fe Volk 2000| |Miniati (2 003)). Here the volume -averaged IC and bremsstrahlung spectra extend 



to GeV-TeV energies and may dominate the 7-ray domain, depending on the strength of the 7r°-decay component. 



3. Hadronic cosmic ray - gas interactions 

Interactions between cosmic ray protons and nucleons of the ICM gas component are very rare and occur on average 
once in a Hubble time in Coma. The number of collisions is usually time-dependent, in particular higher than average 
soon after a (e.g. merger) shock has started to develop, due to newly injected particles. In case of a high cosmic 
ray hadron content 7-rays from the decay of ir° are expected to determine the energy range > 1 GcV. Additionally, 
radiation from the secondary pairs, generated through the decay of charged mesons that are produced by hadronic 
pp-collisions, is expected to contribute to the overall broad band spectrum. The short cooling time scales of those pairs 
radiating in the GHz and hard X-ray to 7-ray band leads to quasi-stationary pair populations at these energies on 
a very short time scale. As a consequence a direct relation between the 7r°-decay 7-ray spectrum and the radiation 
spectrum from the (high energy) secondary pairs is expected. 

The spectral index of Coma's putative relativistic proton distribution nor its normalization has yet been determined 
observationally. Nevertheless some plausible arguments can be found to limit the parameter space. Because cosmic ray 
protons are stored efficiently in galaxy clusters for cosmological times HVolk et al. 19961 |Berezinsky et al. 199T| , the 
radiation from the secondary pairs reflects the injected proton spectrum, and the global proton spectrum should be not 
significantly different from the injected one if uniform injection throughout the cluster is assumed. The structure for- 
mation shock scenario gives injection spectral indices of a p = 2.0 . . . 2.5 for strong shocks (|Mini ati (200 3)|. For merge r 
shocks plunging into the cluster body from the periphery a p can evolve from 2 ... 5 fBerrington & Dermer 2003 ) . 
The normalization of the proton component is limited by three constraints: Firstly, the 7r u -decay 7-rays must not be 
overproduced to violate the EGRET upper limit. Secondly, IC scattering off CMB photons by the secondary pair 
(e ) population produced in pp-interactions leads to a further radiation component that covers the energy range 
from ~ 5 eV (corresponding to e ± of energy ~ 50 MeV) to a few GeV. This component is constrained by the HXR 
flux and EGRET upper limit. The expected non-thermal bremsstrahlung from these secondary e ± lies always below 
the corresponding IC flux level. And lastly, these secondary e ± also emit synchrotron photons, and this leads to a 
constraint imposed by the radio observations. Fig. [2] shows the resulting stationary 7-ray spectra for a p — 2.1,2.3 



and 2.5, calculated using the formalism given in |Pfrommer fe Enfilin (2004) for the 7r°-decay 7-ray production and 



secondary pair production. Wc limited the proton spectrum, assumed to be uniformly injected throughout the cluster, 
to 10 6 GeV since higher energetic protons are difficult to confine within the cluster size l|Colafrancesco &: Blasi 1 998). 
The use of gas and proton density profiles as applicated in e.g. |Blasi (1999)| instead of the volume-averaged param- 
eters leads to only minor changes in the 7r°-decay 7-ray intensity for the here considered volume of Coma (with an 
effective radius of ~ 330 kpc). Above ~ 1 TeV photon absorption due to photon- photon pair production in the cosmic 
infrared-to-optical background radiation field must be taken into account. For this correction we used the background 



models in Aharonian (2001) 



Proton energy densities u p are calculated in the following from the proton spectrum above the threshold for hadronic 
pp-collisions, and are compared to Coma's thermal energy density Mtherm ~ 3.8 ■ 10~ n erg/cm (for kT e — 8.2 keV, a 
thermal e~ density of 10~ 3 cm~ 3 and a 4 He mass fraction of 0.24). The synchrotron flux in the MHz-to-GHz regime 
from the secondary pairs is dependent on u p and a p as well as on B. 

For a p ~ 2.4 and X p = Up/wtherm ~ 20% the radio data are explainable by synchrotron emission from secondary 
in a volume-averaged magnetic field of 0.15/zG if the steepening of the radio spectrum at high frequencies is disregarded, 
in agreement with Blasi k, Colafrancesco (1999)[ Dolag fc Enfilin (2000)| If the steepening of the > 1 GHz radio data is 



taken into account, obviously the synchrotron flux from the secondary pairs must lie below the GHz-radio observations. 
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Fig. 1. Broad band continuum spectrum of Coma. The radio data and the best-fit spectrum at source (corrected for self- 
absorption) are taken from TH03 (not corrected for the thermal Sunyaev-Zeldovich effect). The dotted line represents 
the CMB field corrected for the thermal Sunyaev-Zeldovich effect using a y-parameter of 0.75 • 10 -4 l|Enfilin 2 002). 
The IC and non-thermal bremsstrahhmg fluxes are calculated for field strengths B = 0.1, 0.68, 1.9, 6/1G (from upper 
to lower curves) using an exponential e~ distribution (with adjusted Qq, p c to fit the radio data), and rii — 10~ 3 cm~ 3 . 
They are extended to lower energies assuming the synchrotron spectrum follows a power law down to at least 10 _9 eV. 
The hatched regions in the X-ray domain represent the data from PDS/BeppoSAX I Fusco-Femiano et al. (2004) I, 
flRephaeli et al. 1999D and EUVE l|Lieu et al 



HEXTE/RXTE 
within a ring of 21' and 18' 



1999)) . The RXTE and EUVE fluxes are integrated 
respectively, while the PDS data include only fluxes within 8' from the cluster center. 



In fact, we find that these high frequency radio data place the most stringent constraint on the proton energy content 
in the Coma Cluster. The resulting upper limits for the relativistic hadronic energy density of X p < 3% . . . 0.009%, 
X p < 8%...0.01% and X p < 28%... 0.07% (assuming B = 0.1...2/iG) for a p = 2.1, 2.3 and 2.5, respectively, 
are significantly lower as used in structure formation triggered acceleration scenarios. For example, the model of 



Miniati (2003) required ~ 34% of the thermal energy in form of cosmic ray ions for B = 0.15//G, and ~ 4% for 
B = 0.5/xG with a proton spectrum a p ~ 2 to explain the radio halo emission as originating from the secondary pairs. 
Our cosmic ray limits are also lower than the limits derived from Pfrommer &: Enfilin (2004) (X p < 45% . . . 25% for 
a p = 2.1 . . . 2.5) which solely relied on the EGRET upper limit constraint. For the case a p = 2.1 and B — 0.68/iG we 
find approximative equipartition between particles and fields with X p sa 0.05%. Except for proton spectra harder than 
a p < 2.3 we find in all cases the radiation spectra at > 1 keV from the secondary pairs to lie below the corresponding 
photon spectra from the primaries. This is shown in Fig. 2 for B — 0.1/iG, which simultanously gives the most 
optimistic flux predictions at high energies. Below the hard X-ray band IC from both, primaries and secondary pairs 
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Fig. 2. Same as Fig. 1 but the IC and non-thermal bremsstrahlung fluxes are shown only for a field strength B = 0.1/zG. 
The 7r°-decay 7-ray spectra (most right curves) are calculated for a a p — 2.1 (dashed line), 2.3 (dotted line), 2.5 
(dashcd-dotted line) proton spectrum and the normalization of the particle spectra are adjusted to avoid violating the 
EGRET upper limit as well as the integral fluxes in the HXR and radio domain (see text). The required relativistic 
proton energy densities are 3%, 8% and 28% of the thermal energy content for a p — 2.1, 2.3 and 2.5, respectively. The 
corresponding IC and synchrotron fluxes are shown as dashed/dotted lines. 



determine the shape of the volume-averaged spectrum. Depending on the proton spectral index and overall hadron 
content in Coma, a turnover from primaries' to secondaries' dominated IC below the soft X-ray band may occur. This 
is in agreement with the finding of Bowyer & Berghofer (1998) that the non-thermal halo component detected with 
the EUVE may stem from an additional component of low-energy cosmic ray e~ which we interpret as the secondary 
pairs. Independent hints for a EUV emission of secondary pair origin has just been given by Bowyer et al. (2004) 
who found a striking spatial correlation between the EUVE excess and ROSAT thermal hard X-ray flux based on a 
re-examination of the EUVE data. So far the EUVE excess radiation, if considered to be of non-thermal origin, has 
either been interpreted as IC emission from low energy relic e~ IjSarazin 1999al|Atoyan <k Volk 2000| ) or explained by 
a spectral break between the EUVE and HXR radiating e~ (induced by a certain particle injection scenario), while a 
secondary pair origin had erroneously been ruled out (see Bowyer et al. (2004) for a discussion) . 

In the GLAST energy range non-thermal bremsstrahlung, followed by 7r°-decay 7-rays above ~ 0.1 GeV, will 
dominate, similar to the predictions given in Sarazin (1999b) Only for hard input proton spectra IC radiation from 
the secondary pairs will determine the GeV radiation. No 7-ray emission above ~10 GeV is expected for proton 
injection spectra as steep as a p = 5. 
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4. Detectability with gamma ray instruments 

Advances in the spatial and spectral resolving capabilities of current / future high energy instruments suggest to study 
the Coma cluster continuum emission at energies higher than the hard X-rays. A detection at 7-rays may clarify on the 
spectral extent of the cluster's non-thermal emission, may provide constraints on the acceleration processes realized in 
Coma, and also yields more sensible estimates of a galaxy cluster contribution to the extragalactic 7-ray background 
( |Berrington fe Dermer 20031 |Fujita^t^al^ 2003| ICabici & Blasi 2003a Miuiat i 2002|) . Observationally upper limits are 
currently provided by OSSE iRephaeli et al. 19941 ) and EGRET ijReimer et al. 200311 . INTEGRAL with its moderate 
continuum sensitivity 1 permits a chance to detect Coma as a marginally extended source up to a few 100 keV by ISGRI 
(Gold om et al. 2001jl . PICsIT and SPI, however, will not be able to detect the Coma cluster given realistic observation 
times of < 10 7 sec. The current generation of imaging Cherenkov telescopes (IACTS; Weekes et al. (2002) I, in particular 



at northern locations, will reach the required sensitivity only if significantly more than 50 hrs of observation will be 
accumulated. Even then, the excellent resolving capabilities of IACTs cannot be used to its full advantage due to the 
extended character of Coma's emission where IACTs have a reduced sensitivity. This applies in particular to cases 
where the dominant sub-GeV/TeV-emission component is originating from the outskirts, e.g., due to accretion shocks 
UGabici fc Blasi 2004jl . 

AGILE 2 , expected to have a similar performance as EGRET, might be able to verify the EGRET upper limit. 
It is the Large Area Telescope (LAT) 3 , the main instrument aboard GLAST, that has a realistic chance to finally 
detect Coma in continuum 7-rays. With its significantly better spectral and spatial resolution, and up two orders of 
magnitude improved sensitivity compared to EGRET, the 7r°-decay component will be within reach of LAT. Due to 
the similar spatial extent of Coma C and LAT's point spread function at GeV energies, spatially resolved spectral 
information is difficult to gain. Although being photon-limited, LAT will benefit from its wide field-of-view, that allows 
a steady accumulation of exposure throughout the expected mission life time for any observable object in the sky, 
inclusive the Coma cluster. 



5. Conclusions 

The present work considers the role of the recently confirmed steepening of Coma's radio halo spectrum in the GHz 
band for predicted fluxes in the high energy regime. Indeed, we found that the steepening radio spectrum efficiently 
constrains the amount of hadronic cosmic rays through the radiation channel of the secondary pairs produced in the 
decay chain of the hadronically produced charged mesons. The implied upper limits for the hadronic cosmic ray energy 
density range from 0.01%. . . 28% of the thermal energy density, depending on the magnetic field (B = 0.1/iG . . . 2/iG) 
and proton injection spectral index (a p = 2.1 . . .2.5), and are smaller than those used by other works. This might 
have severe implications for the evolution of galaxy clusters, acting acceleration scenarios in cluster of galaxies and 
the origin (secondary versus primary electron scenario) of Coma's radio halo. 

Below the soft X-ray band we found that a turnover from primaries' to secondaries' dominated IC emission may 
occur, depending on a p and the hadronic cosmic ray content in Coma. This is in agreement with the suggestions of 



Bowyer & Berghofer (1998) that the non-thermal halo component detected with the EUVE may stem from an addi- 



tional population of low-energy cosmic ray electrons which could in this scenario be interpreted as the secondary pair 



component. Independent hints for a EUV emission of secondary pair origin has just been given by Bowyer et al. (2004) 
on the basis of a spatial correlation analysis between the EUVE excess and the ROSAT thermal hard X-ray flux. 

The steepening of the GHz radio spectrum leads to a decline of the IC and bremsstrahlung component of the 7-ray 
spectrum already at 1-10 MeV and a few GeV, respectively, depending on the magnetic field. We have shown that the 
current continuum sensitivity of INTEGRAL'S ISGRI at > a few 100 keV for a 10 6 sec observation is insufficient to 
detect even the most optimistic predicted flux from Coma. The situation is even worse for PICsIT and SPI. 

7r°-decay 7-rays may extend Coma's 7-ray spectrum to TeV energies. However, significant limits to its absolute 
flux are imposed by the radio spectrum (see above). This leads to flux limits that are below the point source minimum 
flux after 50 hrs on-source observations reached by modern generation northern hemisphere Cherenkov telescopes like 
MAGIC and VERITAS. The case is even more hopeless for extended sources. 

All predictions presented here are based on the assumption of power-law proton spectra in the Coma cluster. 
Curved proton spectra may be possible as a result of re-acceleration of the confined cosmic ray hadrons in clusters of 
galaxies l|Gabici fc Blasi 2003b|l . and this might lead to corresponding changes in the predicted limits. 

It will be LAT of the GLAST mission that might finally be able to detect Coma in the 7-ray band if the magnetic 
field and/or Coma's hadronic energy content is favorable. 



1 http:/ /www.rssd.esa.int/Integral/A02/ 

2 http:/ /agile. mi. iasf.cnr.it/Homepage/performances.shtml 

3 http:/ /www-glast. slac.stanford.edu/software/IS/glast_lat_performance. htm 
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